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Gravitational Wave (GW) observations of coalescing compact binaries will be unique 
probes of strong-field, dynamical aspects of relativistic gravity. We present a short review 
of various schemes proposed in the literature to test General Relativity (GR) and alterna- 
tive theories of gravity using inspiral waveforms. Broadly these schemes may be classified 
into two types: model dependent and model independent. In the model dependent cate- 
gory, GW observations are compared against a specific waveform model representative of 
a particular theory or a class of theories like Scalar-Tensor theories. Dynamical Chern- 
Simons theory and Massive graviton theories. Model independent tests are attempts to 
write down a parametrised gravitational waveform where the free parameters take differ- 
ent values for different theories and (at least some of) which can be constrained by GW 
observations. We revisit some of the proposed bounds in the case of downscaled LISA 
configuration (eLISA) and compare them with the original LISA configuration. We also 
compare the expected bounds on alternative theories of gravity from ground-based and 
space-based detectors and find that space-based GW detectors can test GR and other 
theories of gravity with unprecedented accuracies. We then focus on a recent proposal 
to use singular value decomposition (SVD) of the Fisher information matrix to improve 
the accuracies with which Post-Newtonian (PN) theory can be tested. We extend those 
results to the case of space based detector eLISA and discuss its implications. 
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1. Introduction 

General Relativity (GR) is the most studied and well-tested theory of gravity. It has 
been tested to extremely good accuracies by various solar system tests (weak-field 
tests) and binary pulsar observations (strong-field tests) (see Ref.l^^ for reviews). 
Despite these successes, there is motivation to test the theory in regimes which are 
not accessible at present. Firstly, there is no reason why GR is the fundamental 
theory of gravity at all scales. Though solar system and binary pulsar observations 
have tested GR to very good accuracies, there can be deviations from GR in the 
very strong field regime which may not be seen in the solar system or binary pulsar 
observations. There are unresolved problems of singularities in the theory and a 
complete quantum theory of gravity still remains to be a distant goal. Finally, even 
if it is the correct description of gravity, one would like to quantify the statement. 

1.1. Tests of GR and Alternative theories of gravity using GW 
observations 

Gravitational Waves from compact binaries consisting of neutrons stars (NS) and/or 
black holes (BH) are unique probes of strong field aspects of relativistic theories of 
gravity (See Ref.l^ for a review of the field of GW astronomy. See Ref.^^ for an 
overview of electromagnetic counterparts associated with compact binaries). The 
late stages of the evolution of the binary, from which observable GWs are expected, 
are highly relativistic. The velocity of the binary, during the late stages, can be very 
close to the speed of light as opposed to binary pulsars which have speeds ^ 10^'^c. 
Hence GWs from inspiralling compact binaries directly carry the signatures of highly 
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nonlinear, strong-field dynamics of the space time. Hence it is no surprise that GWs, 
once detected, may be used to test some of the fundamental predictions of GR (and 
alternative theories of gravity). 

Evolution of a compact binary system has three phases: adiabatic inspiral, 
merger and ringdown. The inspiral phase is characterized by the adiabatic approx- 
imation ^ ^ 1 where lo is the frequency of the binary and uj the rate of change of 
frequency under GW radiation reaction. This phase can be analytically modelled 
using post-Newtonian (PN) approximation to GR and gravitational waveforms can 
be calculated with very good accuracy (see ReiW for a review of PN theory). The 
highly nonlinear merger phase, where the two compact objects merge to form a sin- 
gle compact object (most likely a BH), can only be modelled through high accuracy 
numerical relativity (see RefEll^ for a review of numerical relativity in the context 
of GWs.). The newly formed compact object loses its asymmetries and settles down 
to an axially symmetric geometry emitting ringdown gravitational waveforms. This 
phase of the binary evolution can be modelled very well using BH perturbation 
theory (see Ref.!^?^. 

The data analysis of coalescing compact binary signals are based on the tech- 
nique of matched filtering (see Ref.*^ for a review of the GW data analysis in the 
case of Gaussian noise.). In matched filtering one cross-correlates the detector out- 
put with a bank of templates (which are copies of pre-calculated gravitational wave- 
forms with different parameters) looking for correlations exceeding certain threshold 
indicating the presence of a signal. As we just saw, all the three phases of binary 
evolution can be accurately modelled analytically or numerically using GR and 
waveforms can be computed. This prior computability of gravitational waveforms is 
central to the detection and parameter estimation of GWs from compact binaries. 

1.2. Classification of the tests 

There have been many proposals in the literature which discuss various types of 
tests that can be performed with GWs. These range from testing some particular 
nonlinear effects of GPSE^IES to very generic tests of alternative theories of gravity 
under a unified frame work.'SHUl wish to (broadly) classify various proposals 
to test GR and alternative theories of gravity into two categories. One where the 
GW data is compared against a particular model or theory of gravity, which we call 
model dependent tests and the second, which are model independent which means 
the GW data is compared against a generic waveform model with free parameters 
which are different for different theories of gravity which can be constrained from 
GW observations. For example, proposals to bound the Brans-Dicke parameter 
using GW observations (e.g. in Ref^HEll) would fall under the model dependent 
test category whereas parametrised tests of PN theorjEHHl and parametrised post- 
Einsteinian frameworkl^^ would belong to the model independent tests category. 
Before we go on and discuss these proposals, it is good to have an idea of the PN 
waveforms, that are used and also about various detector configurations which are 
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already constructed, undergoing upgrades or being proposed. These two aspects are 
discussed in the next two subsections. 



1.3. PN waveforms from inspiralling compact binaries in circular 
orbit 

Using PN approximation to GR, the response of the gravitational wave detector to 
an incident radiation from a source at a luminosity distance is given by 

h{t) = ^[nMF{t)f\osm, (1) 

where Ai = rp^^M is the called the chirp mass of the system (Af being the total 
mass mi + m2 and 77 = "^^'^ the symmetric mass ratio), F(t) = {2TT)~^d^{t)/dt 
is the instantaneous frequency of the radiation, < C < 1 is a dimensionless geo- 
metric factor that depends on position and orientation of the source and its average 
over all orientations is C = 2/5. This waveform is referred to as pattern averaged 
waveform within the restricted waveform approximation, where the PN corrections 
to the phase of the waveform are accounted to the highest accuracy keeping the 
amplitude of the waveform to be leading Newtonian. This may be contrasted with 
the amplitude corrected full waveform where PN corrections to the amplitude are 
also included.'iSHlI! The phasing formula is currently available till 3.5PN {{v/cY) be- 
yond the leading quadrupolar ordeil^^HSIl for nonspinning binaries moving in circular 
orbits. 

Since matched filtering is usually performed in the frequency domain, it is con- 
venient to obtain the Fourier transform of the waveform above, which is efficiently 
done using stationary phase approximation. The frequency domain waveforms read 

^(/)=^r'/'exp[2z^P(/)+*^] , (2) 

with the Fourier amplitude A and phase ^'(/) given by 

C 



A 



*(/) = 2nft, + ^c + Yl [V^fe + In /] (3) 

k 

Here tc and <i>c are the fiducial epoch of merger and the phase of the signal at 
that epoch. There are eight PN coefficients that appear in a 3.5PN phasing formula 
(including the logarithmic terms), which are given byP^^^^ 

A^w^iTTMY'^-^y^a,, (4) 
256 ?7 

i^ki = w^{7rM)('^~''>/^a,i, (5) 
255 rj 

where explicit expressions for and aki can be found, for example, in Eq. (4)-(13) 
of Ref.^^ (which is a slightly modified form of the same given in RefsESESl). 
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Note that the above waveform is specialized to a class of sources which move in 
circular orbits and have no spins. Though the first assumption is somewhat general 
(since gravitational wave radiation reaction is known to circularize the orbit of 
the binary) second assumption of neglect of spins is less general. Waveforms which 
account for the spins and spin induced precession can be found, for example, in 
Ref.^ Discussion of eccentric orbit binaries can be found in Ref.'^^Mzl 

1.4. Various GW detector configurations 

There are different interferometric GW detector configurations, some are already 
operational and undergoing upgrades and some are at the level of proposals whose 
design studies are going on. In the paper we will consider advanced LIGO (aLIGO|^ 
as a representative of the second generation ground-based interferometric GW de- 
tectors. (There are other detectors like advanced Virgo and KAGRA which are 
under upgrade or construction.) We will also consider a proposed 3rd generation 
ground-based interferometer called Einstein Telescope (ET) For space-based de- 
tectors, we will consider the orig inal LISA0 configuration which has recently been 
downscaled to eLISA/NGO.— These are still at the level of design studies and a 
pathfinder mission for LISA is approaching its completion expecting to be launched 
in 2015.'^ We will also consider the proposed space mission DECIGO for which 
design studies are being carried out (see the review by Kent Yagi in this issu^5l!)_ 

The ground-based detectors are sensitive to high frequency GWs in the frequency 
range 1 — 10* Hz, where as LISA/eLISA will be sensitive in the low frequency band 
10~* — 0.1 Hz. DECIGO will have the sensitivity to GWs with frequency between 
10~^ — lOOHz (hence the name Deci Hz Interferometric Gravitational wave Obser- 
vatory) . Together these detectors span a frequency range of 8 orders of magnitude 
between 10"'' - 10* Hz. 

In addition to these detectors there are space-based detectors like ASTROD- 
GWj^hich are currently in the initial stages of its design studies.!^ ASTROD-GW 
will be sensitive to the frequency band between 100 uHz-lmHz. The potential of 
ASTROD-GW to test GR and alternative theories of gravity needs to be explored 
in detail. 

This article is organized in the following way. Sec.|2]reviews some of the model in- 
dependent tests of theories of gravity which includes the scalar-tensor theories, mas- 
sive graviton theories and Dynamical Chern-Simons theory. Proposed generic tests 
of gravity are reviewed in Sec. |3] These include the Parametrised Post-Newtonian 
Tests, Parametrised Post Einsteinian framework and Bayesian methods of perform- 
ing generic tests of strong field gravity. Impact of singular value decomposition for 
these tests are demonstrated in the context of eLISA in Sec. ID Conclusions are 
given in Sec. [5| 



^LIGO is an acronym for Laser Interferometric Gravitational wave Observatory. 
''LISA stands for Laser Interferometer Space Antenna. 

'^ASTROD is an acronym for Astrodynamical Space Test of Relativity using Optical Devices 
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2. Model dependent tests of alternative theories 

In this section, we discuss three important classes of alternative theories: scalar- 
tensor theories, massive graviton theories and dynamical Chern-Simons theory. We 
review literature and summarize the expected constraints on these theories from 
the future ground-based and space-based GW detectors. 



2.1. Scalar-Tensor theories 

Perhaps the simplest and most popular alternative theories of gravity are the scalar- 
tensor theories. In such theories, gravitational interaction is mediated not only by 
the metric tensor but also a scalar field which couples universally to matter. One 
well-studied example of scalar-tensor theories is Brans-Dicke (BD) theory.!^ It will 
be interesting to know what kind of bounds can be obtained on the scalar-tensor 
type theories from GW observations. In the subsections below, we discuss the various 
proposals in the literature discussing the possible bounds on scalar-tensor theories. 
The distinct feature of scalar-tensor theories is that they predict dipolar gravita- 
tional radiation as opposed to GR where the leading order gravitational radiation is 
at quadrupolar order. This would modify the GW phasing formula where a dipolar 
term will appear. One can ask the question how well can GW observations constrain 
this dipolar term. 



2.1.1. Possible bounds on Brans-Dicke theory 

WilP^ argued that the relative importance of the scalar field, in the gravitational 
phasing formula, can be written down in terms of a Brans-Dicke parameter ojbd 
General Relativity is recovered in the limit of large ujbb ('^bd — > oo). The impor- 
tant difference in gravitational wave emission within this theory is the presence of 
dipolar gravitational radiation. WilP^ argued that the PN phasing formula within 
the theory can be written parametrically introducing a (leading order) dipolar term 
in the phasing formula (see Eq. (3.19) o^^^. 

The practical difficulty to obtain a good bound on the BD parameter is that 
gravitational waveforms of BD theory and GR are identical for binary BH systems. 
Further, for binary NS (BNS) systems, the dipolar term is proportional to the square 
of the difference in sensitivity parameter^ which is for a binary BH (which is a 
consequence of the no hair theorem). For binary NSs the sensitivity parameter can 
be between ^ 0.05 — 0.1 and for NS-BH systems ^ 0.3. Thus NS-BH systems are the 
most interesting candidates to bound wbd- Using the matched filtering procedure, 
one can obtain the bounds on this parameter from GW observations of NS-BH 
systems. Fig. [T] compares the projected bounds on ojbd from aLIGO and proposed 



"^This parameter can be a function of the scalar field for generic scalar-tensor theories. 
•^Sensitivity parameter is a measure of self-gravitational binding energy per unit mass (see Sec. 3.3 
of Re^^ 
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ET detectors. The best bounds ~ 10^ would come from a NS-BH systems where 
mass of the BH is less than ^Mq . We also note that if ET can observe binary NSs 
with masses (1.4 + 2)Mq with a sensitivity parameter — 0.1, the projected bound 
on wbd ~ 10'*, if the source is at 100 Mpc. 

LISA may be able to yield interesting constraints on the BD parameter as it 
is sensitive to NS-BH binaries with BH mass ^ 400 — IQ^Mq, if such systems are 
close-by (at a distance of ^ 20 Mpc, which will have an SNR of ^ 10) in which 
case it may yield a bound comparable to or better than Cassini measurements.'^ 
Our calculations, similar to Berti et al.,'^show that bounds expected on ojbd from 
eLISA could be worse by a factor ranging between 3 (for 1 — 4OOAf0) to 6 (for 
(1 — IQ^Mq). On an average the bounds expected from eLISA would be ~ 10*. 

The best bet for bounding BD parameter could be the future probe like DECIGO 
which has sensitivity in the frequency range 10~^ — 10^ Hz which can observe NS- 
BH binaries where the BH is of stellar mass range IQMq). Sensitivity to the low 
frequency GWs as well as the stellar mass range of sources make DECIGO the ideal 
GW detector to yield constraints as high as ~ 10^ as found in Ref.'^ They further 
argue that if they fold in the effect of population of sources DECIGO can observe 
(roughly 0(10*)), the bounds can become as high as 10^. 

Table, [l] summarizes the expected bounds from aLIGO, ET, eLISA, LISA and 
DECIGO. 

Table 1. Projected bounds on ojbd for various detector configurations. These are compared against 
the already existing solar system bounds from the Cassini experiment .ESI Note that an SNR of 10 for 
1 — 400Mq for LISA corresponds to distance of ^ 20 Mpc (and even less for eLISA). For ET, a seismic 
low frequency cut-off of IHz is assumed. 



Detector 


System 


Specification 


Expected bound on wbd 


Reference 


aLIGO 


{lA+b)MQ 


100 Mpc 


~ 100 


Arun irj^ 


ET 


(lA+b)MQ 


100 Mpc 


~ 10^ 


Arun iiflZl 


ET 


(1.4-|-2)Mq (BNS) 


100 Mpc 


~ 10* 


This paper. 


eLISA 


(1.4 + 4OO)Af0 


SNR=10 


~ 10* 


This paper. 


LISA 


(1.4 + 4OO)M0 


SNR=10 


~ 10^ 


Berti et aPl 


DECIGO 


(1.4-1- 1O)M0 


SNR=10 


~ 10^ 


Yagi & TanakiilSl 


Cassini 


Solar system experiment 




10* 


Bertotti et al^ 



Currently, the best bound that exists for the parameter wbd 10^ is from 
Cassini which is a solar system experiment.!^ From the Table it is clear that 
advanced third generation detectors and space based detectors like LISA would be 
the best bets to beat the solar system constraints on the Brans-Dicke parameter. 

Many groups have worked on understanding various features of the waveform 
that may affect the bound on wbd- This includes the effect of LISA noise PSD,I3S1 
effect of nonprecessing spin^MESl and effect of orbital eccentricity.'^ 

Some of the upcoming/proposed solar system experiments are capable of mea- 
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10° h A 




AdvLIG0:(1 .4+5) M at 300 Mpc 



■"^ 5 10 15 20 

BH Mass (Mg) 

Fig. 1. Projected bounds on cjq from the proposed Einstein Telescope for a NS-BH system as 
a function of BH mass (NS mass is assumed to be IAMq). These are compared against the 
best bounds that you expect from aLIGO and also the existing solar system bound from Cassini 
experiment. (Figure is taken from.ES) 

suring the time delay parameter 7 to unprecedented accuracies.'^ For example 
ASTROD-I mission may be able to bound 7 at the level of 3 x 10^^ which is 
three orders of magnitude better than the present bound. This will also lead to very 
accurate constraints on BD parameter, albeit, from non-GW observations. 

2.1.2. Possible hounds on generic scalar-tensor theories 

There have been studies about the bounds that GW observations can put on generic 
scalar-tensor theories. Damour and Esposito-Farese compared the probing power of 
(ground-based) GW observations and binary pulsar observations in the context of a 
two-parameter family of scalar-tensor theories^^ and concluded that binary pulsar 
observations are better probes of strong-field radiative aspects of relativistic gravity. 

More recently, the possible bounds on generic dipolar radiation was studied in 
Ref.,!^ where the bounds on a dipolar phasing parameter /3 (which is related to wbd 
in the case of BD theory) and an amplitude parameter a was obtained for aLIGO 
and ET configurations. A parametrised gravitational waveform in Fourier domain 
(given in Eq. (6) and (7) of that paper) forms the basis of their analysis. Here based 
on the same formalism we extend the results to the case of space-based detectors, 
especially eLISA and LISA detectors. 

Fig. [2] presents the bounds on a and /3 expected from eLISA and LISA config- 
urations. The typical bounds from various detectors are compared in Table [2] The 
table indicates an order of magnitude improvement in the bound on the parame- 
ters using classic LISA configuration, whereas the bounds from ET and eLISA are 
comparable. It should be pointed out that the estimates are obtained using Fisher 
matrix formalism which is valid in the limit of high SNR. 
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Table 2. 
Projected bounds on a 
and /3 parameters describ- 
ing generic dipole GW radi- 
ation.EH For ground-based 
detectors, sources are as- 
sumed to be at 300 Mpc 
and for eLISA and LISA 
sources are assumed to be 
at 3 Gpc. 

a /3 



aLIGO 


10-2 


10-5 


ET 


10-3 


10-'3 


eLISA 


10-3 


10-*^ 


LISA 


10-* 


10-7 



More recently, Berti et al!^ studied possible bounds on the ojBD-like coupling 
parameter if the scalar field has a small mass m^. They showed that compact bi- 
nary observations can bound the combination They obtained (-^|=)(^) < 
10-l^ 10-16, 10-19 eV for aLIGO, ET and eLISA, respectively, where p denotes the 
signal to noise ratio of the observation. 




Mass of the binary (M„) 



Fig. 2. Projected bounds on a and /3, describing generic dipolar gravitational radiation using 
(classic)LISA and eLISA detectors, as a function of the mass of the binary. Sources are assumed to 
be at a distance of 3 Gpc, irrespective of the mass. The errors were obtained using Fisher matrix 
calculation. 



We would like to point out that in our eLISA calculations, we have assumed the 
upper cut-off frequency of eLISA to be O.lHz. The bounds on various scalar tensor 
parameters obtained this way would be worse than those obtained assuming an 
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upper cut-off frequency of IHz. While comparing the eLlSA numbers of this paper 
with that of Yagi)^ it should be kept in mind that Yagi has assumed the upper 
cut-off frequency for eLISA to be IHz and hence the bounds he obtained would be 
better than those reported in this paper. 

2.2. Possible bounds on Dynamical Chern-Simons Gravity 

Chern-Simons (CS) modified gravity theory is one of the most interesting extensions 
of GR where gravitational field is coupled with a scalar field through a parity- 
violating CS parameter (see^for a review). In the dynamical version of the theory, 
the scalar field is treated as a dynamical field. This type of modification naturally 
appears in string theory and loop quantum gravity. Constraints on dynamical CS 
parameter is rather weak because it interacts only with gravity and it reduces to GR 
in the weak-field limit and has the same post-Newtonian parameters as GRUSMHI 
Recently Yagi et al. obtained the leading correction to the GR phasing formula from 
CS type modification.'^ This term occurs at 2PN and they estimated the expected 
bounds on this parameter from GW observations. The CS parameter ('-'^/''scales as 
distance and it tends to GR in the limit ^ ^ 0. Hence the smaller the bound on 
the better. They estimated typical bounds of 0(10— 100)kms with advanced ground- 
based detectors which are better than the expected bound O{10^ — 10®)kms from 
LISA. This is because the expected bounds are of the order of the length scale of the 
target system which in turn is 0{M) where M is the total mass of the system. Since 
LISA is sensitive to super massive black holes, the bounds are weaker. Its worth 
noting that the expected bounds from ground-based and space-based GW detectors 
are 6-7 orders of magnitude stronger than that from solar system experiment, which 
is C^/"* < C'(10«)kms.Eal 

2.3. Massive graviton theories 

Massive graviton (MG) theories are a class of theories where gravitation is propa- 
gated by a massive field as a result of which GWs travel with velocities different 
from the speed of light unlike GR. In this case the velocity of propagation will 
depend up on their frequency as (-^)^ = 1 — (^)^ and the effective Newtonian 

potential may be modelled to be of the Yukawa form cx r~^e where Xg is the 
gravitational compton wavelength. WilP^ argued that, if graviton is massive, it will 
result in different propagation speeds of the low frequency and high frequency GW 
components of an inspiral waveform resulting a modification of the GW phasing pre- 
dicted by GR. This modification can be parametrised in terms of Xg which hence 
can be bounded from GW observations which translates in to bounding the mass 
of the graviton. General Relativity is recovered, obviously, in the limit — >■ oo. 

Many subsequent works improved this model incorporating various features of 
the waveform. Effect of (non-precessing) spins was studied in the context of space 
based LISA detector by Berti et Effect of spin precession was investigated in 
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detail by Stavridis and WilP^ and also by Yagi and TanakapSJ It was found that the 
nonprecessing spins deteriorate the bounds due to large correlation of the massive 
graviton parameter with the spin parameters. But spin induced precession and 
the resulting modulations of the waveform makes up for the deterioration. Effect 
of amplitude corrected PN waveforms was studied by Arun and WilP^ and they 
found an improvement roughly by an order of magnitude, for ground-based and 
space-based detectors for high mass sources (see Fig. 1 of Ref.l^. Our calculation, 
following Ref.,!^ shows that eLISA is expected to yield bounds on Xg which is ~ 10^^ 
kms. The use of FWF can improve the bound by a factor of 7 for high mass systems 
for a mass ratio of 10. Recently Huwyler et al. studied the effect of incorporating 
higher harmonics and spin precession in constraining alternative theories of gravity. 
They found, for precessing case, higher harmonics results in an improvement 1.6 
times compared to the restricted waveforms.!^ 




lu n p 4 fi 

10 10 10 10 

Mass of MBH binary (M^) 

Fig. 3. Comparison of expected bounds for aLIGO, ET and LISA detectors using amplitude cor- 
rected full waveform. Figure taken from RefEB Bounds roughly improve by an order of magnitude 
in going from aLIGO to ET and ET to LISA. 



Effect of orbital eccentricity on the estimation of the massive graviton parameter 
was investigated by Yagi and Tanaka.l^ They found that including (nonprecessing) 
spins and orbital eccentricity results in deterioration of the bound where as invoking 
(simple) precession improves the bound almost by an order of magnitude. Keppel 
and AjitlP^ studied the effect of incorporating the merger and ringdown phases of 
the binary BH waveforms using a phenomenological parametrisation of the numer- 
ical relativity waveforms obtained in Ref.-'* They argued that the inclusion of the 
merger and ringdown parts can lead to drastic improvement for the ground-based 
and space-based detectors. Berti et al. studied^ the effect of population on the 
bounds on Xg and found that that the bounds improve proportional to the square 
root of number of observations. The best expected bounds for different detector con- 
figurations and different waveform models are summarized in Table [3] below. The 
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best bound on Ag would be - lO^^, lO" and lO" kms for aLIGO, FT and LISA, 
respectively. Recently Mirshckari et al!^ proposed a generalization of this idea to 
capture any generic Lorcntz invariance violation which leads to modified dispersion 
relation and discussed the possible constraints from future GW observations. For- 
mulation of the problem from a Bayesian view point was presented in Ref by Del 
Pozzo et al. This method is very handy while implementing the scheme directly in 
the parameter estimation pipelines. 

Table 3. Expected 'best' bounds on \q (in kms ) for various detector configurations and waveform features. FWF refers to 
the amplitude corrected PN waveforrril£ll22l22l2£l (without spins) and IMR is the phenomenological inspiral-merger-ringdown 
waveform obtained by Ajith et al.l^For ground-based detectors aLIGO and ET sources are assumed to be at 300 Mpc and 
for eLISA, LISA and DECIGO, sources are assumed to be at 3 Gpc. Effects of spin, eccentricity etc cannot be read-off from 
the table as we have taken the best possible bounds for each waveform model. 

non-spinning FWF nonpreccssing precessing eccentricity IMR 

aLIGO 10i2(Refgl) ioi3(Ref.Eij) - - - loi3(Refg 

ET 10i3(Ref!S^ 10i4(Ref.ll4 - - - 10i4(Ref.i3^ 

eLISA 10l5(This paper.) 10i5(xhis paper.) - - 

LISA 10i6(Ref!3 loi6(Ref.ll£l) 10i6(RefS 10i''(RefElMl) 10^'' (Ref.^ IQ^T (Rel^ 

DECIGO 1015 (RefB3) - lO^^ {KeiW^ - lO^^ (Ref.EI^ 



3. Model independent Tests 

Having discussed various model-dependent constraints that may be obtained on 
various alternative theories of gravity, this section is devoted to discussing generic, 
model independent tests of GR and alternative theories of gravity. We discuss the 
proposals in chronological order. 

3.1. Parametrised Tests of Post-Newtonian Theory (PTPN) 



As described in Sec. 1.3 3.5PN phasing formula for nonspinning binaries in cir- 
cular orbit is characterized by eight PN coefficients which are predictions of PN 
approximation to GR. An alternative theory would, in principle, predict a different 
set of phasing coefficients. Hence, accurate measurements of these phasing coeffi- 
cients can be interpreted as tests of PN approximation to GR. Thus if one uses 
a parametrised phasing formula, where every term in the PN phasing formula is 
treated as an independent parameter, GW observations can determine these free 
parameters (within the observational error bars) leading to a test of GR.^^^ Since 
eight of these coefficients are fully characterised by two component masses, one can 
perform a consistency test between these PN coefficients in the mass plane of the 
binary (see Fig. 2 of Refl^). 

However, it turns out that there are large degeneracies between these parameters 
which are essentially functions of the two component masses of the binary. We are 
actually using eight parameters to fit a data which is actually characterised by 
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two parameters leading to strong correlations among the parameters l^^l However, 
within these caveats, Ref.^ showed that for a narrow range of masses between 
10^ — lO^Af0 in the LISA band, all of the eight PN parameters can be estimated with 
(1 — cr) relative errors < 1 and the first three PN parameters can be estimated with 
accuracies ~ 10~^, 10~^, 10~^ respectively. As is obvious, though a clean procedure, 
it is not applicable to many sources especially those which ground-based detectors 
can observe. 

This scheme was modified in Ref.l^ where the authors proposed to use only three 
(out of eight) parameters as independent, parametrising the rest of the parameters 
in terms of any of the two PN parameters and perform consistency tests in the 
mass plane of the binarjj^ Advantage of this scheme is that we are dealing with 
only three free parameters which drastically improves the effectiveness of the test. 
The best choice would be to use the first two PN coefficients as basic parameters to 
be used for parametrising all other PN terms except one which may be treated as 
a test parameter. Keeping the basic parameters the same, one can change the test 
parameters yielding different tests (not all are independent). They estimated that 
all the 6 higher order PN parameters, when used a test parameters, can be estimated 
with 1 — (7 relative errors less than unity for a mass range 1 — IOMq for ET-type 
third generation detectors. Almost for the entire mass range, LISA can determine 
them with much better (1-2 orders of magnitude) accuracies, the relative accuracies 
being - 10"^ foj. iq5 _ io6Mo(see Fig. 2 of Ref^. However, with aLIGO not all 
of these test parameters can be estimated with good accuracy. Mishra et al. studied 
the effect of using amplitude corrected PN waveforms on this proposaLl^Sl Employing 
the best available inspiral waveforms with 3.5PN phase and 3PN amplitude, they 
showed that errors on the phasing parameters, especially at high masses, are reduced 
because of the use of amplitude corrected waveforms. While aLIGO can test only 
up to 1.5PN term in the phasing, ET can test the consistency between various PN 
coefficients over the mass range 11 — MMq. Effect of the seismic cut-off frequency 
on these errors and the systematic bias due to spins were also studied. 

Recently Pai and ArurPH showed that starting from the proposal in Ref.,'^ use 
of singular value decomposition can help one obtain the most dominant three new 
parameters which accounts for most of the information. This method of finding the 
three new parameters is more natural. This also brings in aLIGO into the picture for 
which all of these three new parameters can be estimated with a 1 — a relative errors 
10^^ — 10^^ for low mass systems. Accuracies are even more spectacular (^ 10~^) 
for low mass binaries in the ET band. Extension of this proposal to estimate the 
bounds on the new phasing parameters for cLISA and LISA is discussed in Sec. |4] 
below. 



Since two phasing coefficients are necessary to determine the mass, we need at least three phasing 
parameters to perform a consistency test. 
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3.2. Parametrised post-Einsteinian Formalism (PPE) 

Parametrised post-Einsteinian framework"'^^ is developed with the goal of quanti- 
fying the theoretical bias due to the use of general relativistic waveforms in data 
analysis for GWs. The authors wrote down a generic parametrised waveform which 
characterizes the departures from GR through the free parameters which are present 
in the waveform and have different values in different theories. It can be considered 
to be a generalization of FTPN framework to consider more general waveforms 
which may have a different PN structure from GR (in terms of the exponents of 
the frequency). Needless to say, this reduces to PTPN for theories which have a 
PN structure similar to that of GR. They also account for the contribution from 
merger and ringdown phases of the binary evolution. The formalism is, in principle, 
applicable to metric theories of gravity which have same weak-field limit as GR but 
deviate when the gravitational fields are strong. 

The first detailed study to quantify the parameter biases using PPE was reported 
in Ref.^ where they used Bayesian inference and model selection with the simu- 
lated data of advanced LIGO /Virgo network and also for the proposed space-based 
detector LISA. The qualitative conclusions they reached are the following. (1) GW 
observations can constrain various higher order PN coefficients much more tightly 
than binary pulsar observations. (2) Parameter estimates can have significant bias 
if GR templates are used to recover signals which are described by an alternative 
theory of gravity. (3) Detection efficiency of GR templates are also seriously affected 
if the data is characterized by a non-GR theory. 

Further, Chatziioannou et al.l^ put forward a model-independent test of GR 
which permits constraints on the GW polarization contents from GW observations 
of inspiralling compact binaries. Considering three specific alternative theories of 
gravity, viz, BD theory, Rosen's bi-metric theory and Lightman-Lee theory, they 
computed the response function in the case of ground-based GW detectors (for the 
inspiral phase). These helped them to extend the PPE waveform to account for 6 
GW polarization states, the maximum a metric theory of gravity can have. They 
also discussed various data analysis techniques which can lead to the bounds on 
the polarization parameters. PPE and its extensions give a very general frame-work 
for writing down generic waveforms and discussing possible constraints from GW 
observations. 

3.3. Bayesian Inference and Tests of Gravity 

Bayesian approach towards parameter estimation of GW signals is a relatively new 
field (see for example Ref.l^ . If we want to compare two theories of gravity given the 
observed data and all the priors that we have in hand, the first step is to compute 
the posterior probability of the two theories using Bayes theorem. This leads to 
the computation of the important quantity Odds ratio which is nothing but the 
ratio of posterior probabilities of the two theories. Odds ratio is also the product 
of ratios of the priors of the two theories {prior odds) and the Bayes factor. The 
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signal model selection, between GR and an alternative theory, is indicated by the 
corresponding Odds ratio exceeding a threshold value which in turn is based on the 
requirement of a given false alarm rate]^ The estimation of the unknown parameters 
of the theory is carried out by first computation of joint probability distribution 
function of the parameters and then marginalising it over all the parameters except 
the parameter of interest. Computation of these quantities is done via numerical 
simulations involving Monte Carlo techniques such as nested sampling!^ 

In Ref.,'^ Del Pozzo et al. recast the problem of bounding the mass of the 
graviton from the point of view of Bayesian model selection and argued that his 
method can be extended to other alternative theories of gravity. A generic test of 
strong field gravity, within the Bayesian framework, was proposed by Li et alE2IIZil 
In these works they used the proposal of Ref.^^ (described earlier) to test the FN 
coefficients of GR recasting it using Bayesian model selection. They did studies of 
a few cases of deviation from GR, including the cases where phasing has a 1.25FN 
term (which means a a power of frequency which is different from GR). Using 
simulated data, they found that the framework was able to detect such deviations 
even if such deviations are not accommodated in the model waveforms, provided 
the deviations near the most sensitive frequencies of the LIGO band are at least a 
few percent (~ 5 rad at 150Hz.) 

More recently Vallisneri put forward a simple formula, based on Bayesian infer- 
ence, to characterize tests of GR.^^ He argued that, in the limit of large SNR and 
Gaussian noise, the detectability of an alternative theory can be fully characterized 
by a single quantity, viz, the Fitting Factoi''^ between the GR waveform and the 
family of alternative theory of gravity of our interest. Based on his scheme, he found 
that aLIGO type second generation ground-based detectors may be able to measure 
deviations from GR waveforms which are 1 — 10% (corresponding to Fitting factors 
between 0.9-0.99). 

4. Singular Value Decomposition and parametrised tests of 
post-Newtonian theory in the context of eLISA 

Recently Pai and Arun proposed^ the use of truncated singular value decompo- 
sition (SVD) technique to minimise the near-singular nature of the Fisher matrix 
(arising due to strong correlations as seen in Ref.^'* ) in the space of FN phasing co- 
efficients. The truncated SVD approach in the context of parametrised tests of FN 
theory, helps in identifying the most dominant phasing parameters (which are lin- 
ear combinations of the original FN phasing coefficients) which can be determined 
with good accuracy. In Ref. lathis method was described in the context of advanced 
ground-based detectors advanced LIGO and Einstein telescope detectors. In this 



sPalse alarm rate, in the case where GR is compared against an alternative theory of gravity, 
would be the fraction of events for which Odds ratio is above the threshold, but the underlying 
theory is GR. 
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section we extend this method to the case of space-based detectors like eLlSA and 
LISA. 

4.1. Singular Value Decomposition of Fisher Information Matrix 

The basic idea is discussed detail in Sec. 3 of.'^ We present a concise qualitative sum- 
mary of the procedure here not entering too much into the mathematical aspects. 
Fisher information matrix is constructed by the noise weighted inner product of 
tangent vectors of gravitational waveform w.r.t various parameters that are treated 
as independent in the problem. In the presence of Gaussian noise and large SNR, 
square root of the diagonal entries of the inverse of the Fisher matrix (covariance 
matrix) gives 1 — u errors on various parameters.!^ If the parameters describing the 
waveform are strongly correlated, it can lead it ill-conditioned Fisher matrix (see 
Appendix B of Ref.'''* for a discussion). 

Truncated Singular Value Decomposition (SVD) is a technique that is widely 
used in data analysis for reducing the dimensionality of the data by transforming 
correlated variables into a set of uncorrelated ones and understanding which of the 
newly constructed parameters are the most significant to describe the data and re- 
express the data in terms of these most significant ones. In our case, as we saw in 
Sec. |3.1[ there is strong correlation among the eight FN phasing parameters when 
all of the are treated as independent. We can apply the scheme of Truncated SVD to 
obtain the most significant parameters which are linear combination of the original 
FN coefficients. Below is a more technical details of this procedure. 

We write the Fisher information matrix in the 8 dimensional FN phase parameter 
space, namely 9 = {ipo, V'2, V'3 j V'5 j V'si ; V'e j V'e; , V'y}: in terms of 0" with a = 1 ... 8 
in the discrete frequency domain as 



F,j - m 



.k 



A/. (6) 



where at fc-th frequency bin /j,, h i{f^.) — dh{fk)/d9i denotes the i ^th component 
of the tangent vector of the GW signal in the intrinsic FN phase space and A/ is the 
width of the frequency bin. We stack all the tangent vectors in a matrix form with 
column index representing the frequency index as H = {hik = ^^^'^'H==^}sxn- 

V S(fk) 

Then the Fisher matrix takes a compact form in terms of H as 

r = 2HH^ (7) 

where the dagger (f ) denotes transposed conjugate of the matrix. 

Though not very surprising result, this compact form of Fisher enables us to 
treat the ill-conditioning of the Fisher which arises due to the correlations between 
the parameters which was one of the main hurdles in the previous work.lSl We note 
that H acts like a building block for the Fisher. It carries the information about the 
spectral content of the GW signal variation in the interferometer's noise in 9 space. 
The natural way to remove the correlations would be to look for those directions 
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in the parameter space where the signal variation is maximum. Technically, these 
are the principle directions or principal components of the Fisher matrix. The SVD 
provides the natural way to arrive at the principal components. 

The SVD of H is H = UnJ^nVH^ where UHmxm a-nd Vn^xn are left and right 
unitary square matrices and flu arranges the singular values in descending order as 
diagonal entries. The singular values govern the contribution of the singular vector 
in a particular direction. Substituting the SVD of H in F, we obtain the SVD of F 
as 

r = UHSu;j, (8) 

where S = 2JlHf^H- 

4.2. Truncated SVD, new phase parameters and relative errors 

For a near-singular matrix, some of the singular values are very small compared to 
the most dominant one i.e. very little information lies along those directions. In the 
SVD of H, the signal variation along those directions are very small or negligible 
which makes the matrix ill-conditioned. This ill-conditionedness in turn gives rise 
to unreliable estimation of errors. We propose to truncate the singular matrix of H 
to the dominant singular values. If r is the number of dominant values, the Fisher 
matrix in terms of truncated SVD is 

r = u^jS*u^ ^ u^ru^j = s* (9) 

where the index t corresponds to the truncation. The U^g^^, gives the singular- 
vectors corresponding to the dominant singular values. The the singular values 
S* is a diagonal matrix with singular entries Sfcfe = 2(ri^)^j, for A; = 1, ... , r. 
The truncation criterion is e < (S^H)kk/i^H)ii with arbitrary e. The maximum 
value of k which satisfies the above mentioned criterion gives the rank r. 

As shown in Eq. Q, the singular vector not only gives the singular values 
of the Fisher matrix but also gives transformation into the new phase parameters 
through linearly combination of 6°'. The new phase parameters are 

0' = V^^e. (10) 

where 9' = {V'l, V's! ■ ■ • : V'rl is a r < 8 dimensional phase vector. The Fisher 
matrix and hence the covariance matrix are diagonal in 6' space. Further, the new 
phase parameters 6' are statistically independent variables with zero mean, variance 
(T^ ~ (S^Q,)"^ which provides the absolute error. 

4.3. Relative errors for eLISA noise configuration 

We apply the above technique for the eLISA and LISA configurations using the 
noise curve given in Eqs (l)-(5) of Refl^aland Eqs (2.28)-(2.31) of Ref.^Sl We ignore 
the orbital motion of the space craft and use pattern averaged waveforms. We retain 
the three most significant new parameters which corresponds to e = 10^^. 
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eLISA Classic LISA 




io' io' io' 10° io' 10° 10* io' 10* io' 

Total Mass In Solar Mass Units Total Mass In Solar Mass Units 



Fig. 4. Errors in the most dominant new parameters (which are linear combinations of the old 
parameters) obtained from SVD in the LISA and eLISA bands as functions of the total mass of the 
binary for equal mass systems rj = 0.25. Distance to the source is assumed to be 3 Gpc irrespective 
of the total mass. Spins of the binary systems are ignored. 

In Fig. |4j we plot the relative errors of the new parameters with the total mass 
for the dominant parameters for eLISA and classic LISA cases. It is obvious that 
cLISA has poorer estimation of the new parameters roughly an order of magnitude 
worse than classic LISA for all the three parameters. This may be attributed to the 
reduced sensitivity of eLISA, compared to LISA, due to change in the configuration. 
Except for the spiky segment, all the errors are less than lO"'^ for eLISA and 10"''^ for 
LISA highlighting the unprecedented accuracies with which these new parameters 
can be extracted by space-based detectors which are sensitive to the low frequency 
GWs. Implications of this for tests of GR are currently under investigation. We give 
below a Table of accuracies with which three leading parameters may be extracted 
by aLIGO, FT, eLISA and LISA. 



Table 4. Typical relative errors -jp^ 

k 

of the new phasing coefficients ob- 
tained via SVD, for k = 1,2,3 for 
aLIGO, ET, eLISA and LISA configu- 
rations. 
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5. Conclusion 

We have presented an overview of various proposed methods to test GR and alter- 
native theories of gravity using GW observations of inspiralling compact binaries. 
We compared, in each case, the ability of various existing, ongoing and proposed 
ground-based detectors and the proposed space-based detectors to constrain theories 
of gravity. It is very clear that one needs to go the low-frequency GW observations 
of space-based detectors to test GR with very good accuracy. Hence missions like 
eLISA/NGO, LISA, DECIGO are extremely important for fundamental physics. 

As a final note, we would like to point out that this is not a comprehensive 
review of the field. Also there are various very interesting proposals using merger 
and ringdown waveforms of the compact binaries which is not discussed in this 
article. 



6. Acknowledgements 

KGA would like to thank the organizers of the ASTROD symposium in Banga- 
lore. We thank Bala Iyer, Wei-Tou Ni and Kent Yagi for useful comments on the 
manuscript. The work is supported in part by AP's MPG-DST Max-Planck India 
Partner Group Grant. 



References 

1. C. M. Will, Theory and experiments in gravitational physics (Cambridge University 
Press, New York, USA, 1981). 



2. C. M. Will, Living Rev. Rel. 9 (2005) 3, |arXiv : gr-qc/0510072 

3. I. H. Stairs, Living Rev. Rel. 6 (2003) 5, |arXiv : astro-ph/0307536 [astro-ph] [ 



arXiv: 0903. 0338 



4. B. Sathyaprakash and B. Schutz, Living Rev. Rel. 12 (2009) 
[gr-qc] 

5. A. Kamble and D. Kaplan, Int. J. Mod. Phys. D 22 (2013) 1341012. 

6. L. Blanchet, Living Rev. Rel. 9 (2006) 4,^gr-qc /0202016 

7. F. Pretorius (2007) |arXiv: 0710. 1338 [gr-qc] I Relativistic Objects in Compact Bi- 
naries: From Birth to Coalescence Editor: Colpi et al. 

8. J. Centrella, J. G. Baker, B. J. Kelly and J. R. van Meter, Rev.Mod.Phys. 82 (2010) 
3069, arXiv: 10 10. 5260 [gr-qc] | 

9. M. Sasaki and H. Tagoshi, Living Rev. Rel. 6 (2003) 6, [arXiv : gr-qc/0306120' 

10. E. Berti, V. Cardoso and A. O. Starinets, Class. Quant.Grav. 26 (2009) 163001, 
jarXiv : 0905 . 2975 [gr-qc] 

11. P. Jaranowski and A. Krolak, Living Reviews in Relativity 15 (2012). 

12. L. Blanchet and B. S. Sathyaprakash, Class. Quantum Grav. 11 (1994) 2807. 

13. L. Blanchet and B. S. Sathyaprakash, Phys. Rev. Lett. 74 (1995) 1067. 

14. K. G. Arun, B. R. Iyer, M. S. S. Qusailah and B. S. Sathyaprakash, Class. Quantum 
Grav. 23 (2006) L37, gr-qc/0604018 

15. K. G. Arun, B. R. Iyer, M. S. S. Qusailah and B. S. Sathyaprakash, Phys. Rev. D 74 
(2006) 024006, lgr-qc/0604067 



16. N. Yunes and F. Pretorius, Phys. Rev. D 80 (2009) 122003, jarXiv: 0909 .3328 
I [gr-qc] I 

17. C. M. Will, Phys. Rev. D 50 (1994) 6058, |arXiv : gr-qc/9406022l 



February 12, 2013 2:37 WSPC/INSTRUCTION FILE ASTROD' article 



20 Arun and Pai 

18. A. Krolak, K. Kokkotas and G. Schafer, Phys. Rev. D 52 (1995) 2089. 

19. K. G. Arun, L. Blanchet, B. R. Iyer and M. S. S. Qusailah, Class. Quantum Grav. 21 
(2004) 3771, gr-qc/0404185, Erratum-ibid. 22, 3115 (2005). 

20. L. Blanchet, G. Faye, B. R. Iyer and S. Sinha, Class. Quantum. Grav. 25 (2008) 
165003, [arXiv: 0802 .1249 [gr^qc]] 

21. A. M. Sintes and A. Vecchio, Detection of gravitational waves from inspiraling com- 
pact binaries using non-restricted post-newtonian approximations, in Rencontres de 
Moriond: Gravitational waves and experimental gravity, ed. J. Dumarchez (Frontires, 
Paris, 2000). gr-qc/0005058 

22. A. M. Sintes and A. Vecchio, Lisa observations of massive black hole binaries us- 
ing post-newtonian waveforms, in Third Amaldi conference on Gravitational Waves, 
ed. S. Meshkov (American Institute of Physics Conference Series, 2000), p. 403. 
[gr-qc /0005059 

23. C. Van Den Broeck and A. Sengupta, Class. Quantum Grav. 24 (2007) 155, 
|gr-qc/0607092 

24. K.^G. Arun, A. Buonanno, G. Faye and E. Ochsner, Phys. Rev. D 79 (2009) 104023, 
farXiv : 0810 . 5336 [gr-qc] 

25. L. Blanchet, T. Damour, B. R. Iyer, G. M. Will and A. G. Wiseman, Phys. Rev. Lett. 
74 (1995) 3515, gr-qc/9501027 

26. L. Blanchet, G. Faye, B. R. Iyer and B. Joguet, Phys. Rev. D 65 (2002) 061501(R), 
gr-qc/0105099, Erratum-ibid 71, 129902(E) (2005). 

27. L. Blanchet, T. Damour, G. Esposito-Farese and B. R. Iyer, Phys. Rev. Lett. 93 (2004) 
091101, gr-qc/0406012 

28. T. Damour, B. R. Iyer and B. S. Sathyaprakash, Phys. Rev. D 63 (2001) 044023, 
|gr-qc/0010009, Erratum-ibid. D 72 (2005) 029902. 

29. T. Damour, B. R. Iyer and B. S. Sathyaprakash, Phys. Rev. D 66 (2002) 027502, 
rgr-qc/0207021 , Erratum-ibid 66, 027502 (2002). 

30. K. G. Arun, B. R. Iyer, B. S. Sathyaprakash and P. A. Sundararajan, Phys. Rev. D 
71 (2005) 084008, gr-qc/0411146, Erratum-ibid. D 72,069903 (2005). 

31. A. Pai and K. Arun (2012) arXiv: 1207. 1943 [gr-qc] 

32. C. K. Mishra, K. G. Arun, B. R. Iyer and B. S. Sathyaprakash, Phys. Rev. D 82 
(2010) 064010, arXiv: 1005.0304 [gr-qc] 

33. A. Gopakumar and B. R. Iyer, Phys. Rev. D 56 (1997) 7708. 

34. T. Damour, A. Gopakumar and B. R. Iyer, Phys. Rev. D 70 (2004) 064028, 
gr-qc/0404128 

35. K. G. Arun, L. Blanchet, B. R. Iyer and M. S. S. Qusailah, Phys. Rev. D 77 (2008) 
064035, , 0711 .0302 , 

36. K. G. Arun, L. Blanchet, B. R. Iyer and M. S. S. Qusailah, Phys. Rev. D 77 (2008) 
064034, 071170250 

37. K. G. Arun, L. Blanchet, B. R. Iyer and S. Sinha, Phys. Rev. D 80 (2009) 124018, 
|arXiv : 0908 ■ 3854 [gr-qc] | 

38. http://www.et-gw.eu/. 

39. P. Amaro-Seoane, S. Aoudia, S. Babak, P. Binetruy, E. Berti et al, Class. Quant. Grav. 
29 (2012) 124016, arXiv: 1202 . 0839 [gr-qc] 

40. P. W. McNamara, Int. J. Mod. Phys. D 22 (2012) 1341001. 

41. K. Yagi, Int. J. Mod. Phys. D 22 (2013) 1341013^ 

42. W.-T. Ni, J. Adv. Space Res. (April 2012) arXiv : 110 4 . 5049 [astro-ph . CD]"] 

43. C. Brans and R. H. Dicke, Phys. Rev. 124 (1961) 925. 

44. E. Berti, A. Buonanno and C. M. WiU, Phys. Rev. D 71 (April 2005) 084025, 
gr-qc/0411129. . 



February 12, 2013 2:37 WSPC/INSTRUCTION FILE ASTROD' article 



Tests of General Relativity and Alternative theories of gravity using GW observations 21 



45. K. Yagi and T. Tanaka, Prog. Theor. Phys. 123 (2010) 1069, arXiv: 0908 .3283 

46. B. Bertotti, L. less and P. Tortora, Nature 425 (2003) 374. 

47. M. Abernathy et al, Einstem gravitational wave Telescope: Conceptual Design Study 
(Document number ET-0106A-10) . 

48. C. M. Will and N. Yunes, Class. Quantum Crav. 21 (2004) 4367, gr-qc/0403100' 

49. K. Yagi and T. Tanaka, Phys. Rev. D D81 (2010) 064008, arXiv : 0906 . 4269 [gr-qcTl 
Erratum tbidD 81, 109902 (2010). 

50. C. Braxmaier, H. Dittus, B. Foulon, E. Goklu, C. Grimani et al, Exper.Astron. 34 
(2012) 181, |arXiv: 1104.0060 [gr-qc][ 



51. T. Damour and G. Esposito-Farese, Phys. Rev. D 58 (1998) 042001, gr-qc/980303lj 

52. K. G. Arun, Class. Quant. Crav. 29 (2012) 075011, arXiv: 1202.5911 [gr-qc] 

53. E. Berti, L. Gualtieri, M. Horbatsch and J. Alsing (2012) arXiv : 1204 . 4340 [gr-qc]] 

54. S. Alexander and N. Yunes, Phys. Rept. 480 (2009) l, ,arXiv: 0907 .2562 [hep-thlf 

55. S. Alexander and N. Yunes, Phys. Rev. Lett. 99 (2007) 241101, |arXiv : hep-th/070 3:565l 
[hep-th] 



56. S. Alexander and N. Yunes, Phys. Rev. D 75 (2007) 124022, ,arXiv : 0704 . 0299 



I [hep-th] I 



57. K. Yagi, N. Yunes and T. Tanaka (2012) arXiv: 1208.5102 [g r-qc] 

58. Y. Ali-Haimoud and Y. Chen, Phys. Rev. D 84 (2011) 124033. 

59. C. M. Will, Phys. Rev. D 57 (1998) 2061, gr-qc/9709011 , 



60. A. Stavridis and C. M. Will, Phys. Rev. D 80 (2009) 044002, a rXiv : 0906 ■ 3602 

61. K. G. Arun and C. M. Will, Class. Quant. Crav. 26 (2009) 155002, 



arXiv: 0904. 1190 



62. C. Huwyler, A. Klein and P. Jetzer (2011) arXiv: 1108. 1826 [gr-qc]_ 

63. D. Keppel and P. Ajith, Phys. Rev. D 82 (2010) 122001, arXiv : 1004 .0284 [gr-qc]^| 

64. P. Ajith, S. Babak, Chen et al, Phys. Rev. D 77 (May 2008) 104017, 0710.2335 

65. E. Berti, J. Gair and A. Sesana, Phys. Rev. D 84 (2011) 101501, arXiv: 1107 .35281 
I [gr-qc]] 

66. S. Mirshekari, N. Yunes and C. M. WiU, Phys. Rev. D 85 (2012) 024041, 
[arXiv: 1110.2720 [gr-qc] 

67. W. Del Pozzo, J. Veitch and A. Vecchio, Phys. Rev. D83 (2011) 082002, 
larXiv: 110 1 .1391 [gr-qc] | 

68. L. Blanchet, B. R. Iyer, C. M. Will and A. G. Wiseman, Class. Quantum Crav. 13 
(1996) 575, gr-qc/9602024 

69. C. Van Den Broeck and A. S. Sengupta, Class. Quantum Crav. 24 (2007) 1089, 
[gr-qc /0610126 

70. N. Cornish, L. Sampson, N. Yunes and F. Pretorius, Phys. Rev. D 84 (2011) 062003, 
|arXiv:1105T2 "088 [gr-qc] 

71. K. Chatziioannou, N. Yunes and N. Cornish (2012) arXiv: 1204. 2585 [gr-qc] 

72. J. Veitch and A. Vecchio, Phys. Rev. D 81 (2010) 062003, a rXiv: 0911 . 3820] 
I [astro-ph . COT| 

73. T. G. F. Li et al, Phys. Rev. D 85 (2012) 082003, arXiv: 1111 . 5274 [gr-qc] I 

74. T. Li, W. Del Pozzo, S. Vitale, C. Van Den Broeck, M. Agathos et al, J. Phys. Conf.Ser. 
363 (2012) 012028, arXiv : 1111 . 5274 [gr-qc] 

75. M. Vallisneri, Phys. Rev. D 86 (2012) 082001, ar Xiv: 1207.4759 [gr-qc]] 

76. T. A. Apostolatos, Phys. Rev. D 52 (Jul 1995) 605. 

77. C. Cutler and E. Flanagan, Phys. Rev. D 49 (1994) 2658. 



